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I Abstract. We study the effects of peculiar velocities on statis- 
^ tical measures of galaxy clustering. These effects occur when 
^^ distances to the galaxies are estimated from their redshifts. It is 
'/"J assumed that the clustering pattern results from the gravitational 
^~~^ instability of initially Gaussian, small-amplitude perturbations 
of a Friedman-Lemaitre cosmological model. Explicit expres- 
Q/ sions are given for an arbitrary density parameter Q of the 
21 model, both when the cosmological constant, A, is zero, and 
l> when the model is spatially flat, il + A/3H^ = 1. 

^^ Kaiser (1987) had analyzed the redshift distortion of the 

^" two-point correlation function. This function determines the 

■!i^ variance of the density field distribution function and can be 

F^ computed using linear perturbation theory. We show here how to 

♦recompute higher order moments in redshift space, paying special 

Q attention to the skewness, or third moment of the density field, 

jH and its Fourier space counterpart, the bispectrum. This calls for 

^ a (weakly) non-linear analysis. 

We rely on a perturbative expansion of particle trajectories 
in Lagrangian coordinates, using the formalism introduced by 
Moutarde et al. (1991) and further developed by Bouchet et 
al. (1992, 1994). This formalism extends to higher orders the 
Zel'dovich first order (i.e. linear) solution (1970). The lowest 
non-vanishing contribution to the skewness comes from the 
first and second-order terms in perturbation theory. Therefore, 
using Zel'dovich approximation would not be self-consistent 
and would yield inaccurate results. We show that a physically 
consistent and quantitatively accurate analysis of the growth 
skewness in redshift space can be obtained from second-order 
Lagrangian theory. 

With practical applications to redshift surveys in mind, we 
also study the effects of spatial smoothing ofthe evolved density 
field. The necessary formalism was developed by Juszkiewicz 
and Bouchet (1991) and Juszkiewicz et al. (1993a). Here we 
give the first complete account of these calculations; we also 
extend the formalism by explicitly taking redshift distortions 
into account. We give analytic expressions for the gravitation- 



ally induced skewness as a function of the power spectrum and 
of il, for a spherical top-hat and a Gaussian smoothing filter. 

We compare our analytical predictions with measurements 
performed in numerical simulations, and find good agreement. 
These results should then prove useful in analyzing large scale 
redshift surveys. In particular, our results, in conjunction with 
the recent suggestion of Fry (1994), may solve a well known 
problem which always arises in conventional dynamical deter- 
minations of the mean density of the universe. Such studies 
produce estimates of il which are coupled with the parameters 
describing the bias in the galaxy distribution. As a result, a 
biased il = 1 model is dynamically indistinguishable from an 
open, unbiased, one. For the first time, it may become possible 
to break this degeneracy, and decouple the estimates of linear 
and non-linear bias from the estimates of il and A. 
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1. Introduction 

The appearance of structures is distorted in redshift space by 
peculiar velocities. At "small" scales, this leads to the "finger of 
god" effect: pattern are elongated along the line-of-sight due to 
their internal velocity dispersion. This is an intrinsically non- 
linear effect. At "large" scales, the effect is reversed: coherent 
inflows lead to a density contrast increase parallel to the line- 
of-sight. Indeed, foreground galaxies appear further than they 
are, while those in the back look closer, both being apparently 
nearer to the accreting structure. 

Bean et al. (1983) and Davis and Peebles (1983) took ad- 
vantage of the small scale distortion of clustering to estimate 
the pairwise velocity distortion. They were followed my many 
others. On the other hand, Sargent and Turner (1977) showed 
that the large scale distortion could be used to measure the cos- 
mological parameter il; Kaiser (1987) pointed out by using a 
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linear analysis that, in redshift space, each wave of wavevec- 
tor k and amplitude S]. is amplified by a factor (1 + /i(il)//|), 
where /i ~ il" *" describes the linear growth of perturbations 
[Peebles 1976, 1980 §14, and Eq. (16) below], and //^ is the 
cosine of the angle between k and the line-of-sight. The true 
power spectrum, P(k) = {\6k\^) in redshift space is replaced 



by 

PAk)=[i + fi(n)filfp(ky, 



(1) 



its Fourier transform, the two-point correlation function ^2, is 
strongly anisotropic in redshift space. This compression along 
the line-of-sight was computed by Lilje and Efstathiou (1989) 
and McGill (1990). Practical improvements to measure il were 
suggested by Hamilton (1992a), and used (Hamilton 1992b) 
to determine il from the IRAS 2 Jy redshift survey (Strauss 
et al. 1990, 1992). Recently, Fisher et al. (1994) measured the 
redshift space ^2 in the IRAS 1.2 Jy redshift survey. Many 
authors preferred though to measure the redshift space power 
spectrum -P^(fc) (see e.g. Feldman, Kaiser, and Peacock 1994; 
Cole, Fisher, and Weinberg 1994, and references therein). 

Most of the analytical calculations performed in this field 
were done using a purely Eulerian approach. Still, Bouchet et 
al. ( 1 992, hereafter B JCP), building on the work of Moutarde et 
al. (1991), developed a systematic perturbative approach from 
a Lagrangian point of view. This work was recently extended by 
Bouchet et al. (1994, hereafter BCHJ), complementing a paral- 
lel effort by Buchert (1994, and references therein). As we shall 
see, it is straightforward to use the Lagrangian approach to in- 
vestigate the redshift distortion effect in the weakly non-linear 
regime. For example, Bouchet et al. (1992) indicated that the 
skewness factor 53 (see below) of the unsmoothed density field 
is nearly invariant under the transformation to redshift space. 
The corresponding demonstration, which uses an "infinitely re- 
mote observer" approximation is given in Bouchet et al. (1994, 
hereafter BCHJ). In this paper we consider a slightly different 
approximation, valid in the limit of large samples, similar to the 
limit considered by Kaiser (1987). It gives a result, nearly iden- 
tical to the one of BCHJ. We generalize the first order power 
spectrum results of Kaiser to second order quantities, such as 
the bispectrum ((5^,(5^,(5^3) J. +k +k =0 ^'^'^ i*-^ '^°™^'" '-'^'^'^'■^'"P^'-' 
53. 

The skewness factor, S3, is the ratio of the third moment 
of the probability distribution function (PDF) of the density 
contrast S = p/ (p) — Ito the square of its second moment (its 
variance) 



S3 



{s^y 



(2) 



Throughout this paper, the symbols (...) denote ensemble av- 
erages, and the redshift space counterpart of a quantity A is 
denoted by A^ . In (Eulerian) linear theory, a normal PDF re- 
mains normal, so that S3 (as well as higher order factors Sn) 
is zero for an initially Gaussian field. Of course, linear theory 
(Eulerian or Lagrangian) is inadequate to describe indicators 
which are intrinsically of higher order. Under the influence 



of weakly-nonlinear gravitational instability, mode-mode cou- 
plings occur and establish a non-zero value for S3, which re- 
mains constant until strongly non-linear effects come into play. 
For any smoothing scale, the ratio S3 describes the asymmetry 
between the positive and negative tail of the PDF. As was shown 
recently by Juszkiewicz et al. (1993b), the low order moments 
of the PDF are enough to describe the weakly non-linear evolu- 
tion of the PDF, thanks to an Edgeworth expansion in powers of 
the variance. The contribution of the third order moment i.e. of 
the skewness itself, can indeed already give a good idea of the 
overall shape of the PDF for ((5^) ^ 0.5 (more precisely, the 

— 1 /2 

tails are reliably described within ((5^) standard deviations 
from (5 = 0). 

Juszkiewicz and Bouchet (1991, hereafter JB) and 
Juszkiewicz et al. (1993a, hereafter JBC) considered the effect 
of smoothing on S3, when the evolved density contrast field S 
is spatially averaged with a weight W, called "window func- 
tion". In the unsmoothed limit obtained by Peebles (1980) for 
an il = 1 model, S3 depends neither on scale nor on the power 
spectrum, and equals to 34/7. But the skewness factor of the 
smoothed density field depends on both the shape of the win- 
dow function and the local slope of the power spectrum. If P(k) 
is not scale-invariant (as in CDM for instance), this induces a 
further scale dependence. For a top-hat window (a sphere with 
sharp boundaries) and a scale-invariant power spectrum with an 
index n in the range — 3 < n < 1, they obtained 



53 = f4(Q-/« 



l)-(rn-3) 



(3) 



and the above result was derived for real space. JBC considered 
Gaussian filters as well. Bernardeau (1993, 1994) generalized 
the result above for higher order S'at , # > 4, in the case of a top- 
hat smoothing filter, while the kurtosis in the Gaussian smooth- 
ing filter case was recently obtained by Lokas et al. (1994). 
These analytical results are fully confirmed by numerical sim- 
ulations (e.g. Bouchet and Hernquist 1992, Weinberg and Cole 
1992, JBC, Juszkiewicz etal. 1993b, Lucchinetal. 1994, Lokas 
et al. 1994). In this paper, we give the first complete account of 
the technique introduced by JB and JBC to estimate the skew- 
ness of the smoothed density field. These calculations, which 
here appear in the Appendix B, were earlier presented by JB 
and JBC in a rather terse form' When necessary, the formalism, 
presented here is also extended to combine smoothing effects 
with the effects of mapping from real space to redshift space. 
This enables comparisons between analytical predictions for 
an initially Gaussian field and practical measurements in the 
observed three-dimensional galaxy catalogs, such as those re- 
cently performed by Bouchet et al. (1991, 1993) in the IRAS 
1.2Jy galaxy redshift catalog and by Gaztafiaga (1992) in the 
CfA and SSRS redshift catalogues. 



Although we never used the the usual hypocritical "it is easy to 
show" phrase, our Letter to the Ap.J. (JBC) may have suggested that 
implicitly. Since it was not so easy at all (in fact it has taken several 
months to derive these results), and since we decided that it is awkward 
to keep distributing hand-written notes to colleagues, asking for details 
of our calculations, we decided to include them in this paper. 
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This paper is organized as follows. In Sect. 2, we recall 
the needed Lagrangian theory results and show how to derive 
the effect of peculiar velocities. We analyze the bispectrum in 
Sect. 3. Smoothing is considered in Sect. 4, where we compare 
our analytical results to numerical simulation measurements. 
Section 5 is the conclusion. 

2. Redshift Distortion in Lagrangian Theory 

We start by recalling some results of BJCP and BCHJ. The 
following section is devoted to the specific case of the redshift 
distortion. 

2.1. Basics of the Perturbative Lagrangian Approach 

We consider a non-relativisticpressureless fluid, in a Friedman- 
Lemaitre background model in the the matter era. During that 
epoch, and at the scales of interest ( ^ lOOMpc), the perturba- 
tions can be treated as a ideal pressureless fluid in the Newtonian 
approximation (cf. Peebles, 1980). The scale factor of the met- 
ric is denoted by a. Throughout, we use comoving coordinates 
X = r /a, where r stands for physical (expanding) coordinates. 
In the Eulerian approach, the primary quantity of the analysis is 
usually the density contrast field b. In the Lagrangian picture, 
each fluid element is tagged by its unperturbed Lagrangian 
(comoving) coordinate q, and the corresponding value of the 
displacement field *. This vector field connects at any time the 
element's Eulerian and Lagrangian position through 



with Vg = {d/dqi, d/dq2, d/dq^). The assumption of irrota- 



tionnal flow leads to the constraints 



V„ 



* 



(1) 



0, 



v„ 



* 



(2) 



:0. 



(10) 



The displacement field ^'^'^ may be related to the initial Eu- 
lerian density contrast Si{x) = 0{s) through the equation (5) 
expanded to linear order 



eV,.^''>(t,q) = 6[x(q)] + 0(e'), 



(11) 



and evaluated at t = ti on the l.h.s and at a; = q on the r.h.s.. 

The n-th order growth rates, </„ depend on the matter content 
of the model. For a background model with A = 0, the fastest 
linear growing modes is identical to the well-known linear 
Eulerian growth rate D(t) of the density contrast (e.g. Peebles 
1980, BJCP), as can be seen from (11). The second order growth 
rate behaves according to 



,2^-^fi-2/«,^ A = 0. 



(12) 



This approximation is better than 0.4% for 0.05 ^ il ^ 3 
(BJCP). We will also consider a flat Universe with A 7^ 0, 
that is il H- A/3H^ = 1. The second order growth factor is then 
given by (BCHJ) 



92 



:fi 



■'''''9I 



n + A/3H^= 1. 



(13) 



x = q + ^(t,q). 



(4) Details of the derivations of these results may be found in BCHJ. 



The Jacobian of the transformation from q to a; permits 
to express the requirement of mass conservation simply as 

p(x)\dx/dq\ d^q = p{q)d^q, or 



S[x{q)] = 



dx 



dq 



(5) 



In the weakly non-linear regime, solutions of the motion 
and field (Poisson) equations can be obtained by means of a 
perturbative expansion of the displacement field 



2,T,(2)^c3,T,(3). 



* = £*^'^ + e^*^^-' + e'-^'^"' + 



(6) 



where e is a formal parameter allowing to conveniently keep 
track of the various orders in the calculations. In the following, 
we need only the first two orders, which are separable: 



>(i) 



.T.(l)/ 



^'"{t,q) = g,{m-^'{q) 



>(2) 



,T.(2)/ 



*^''(i,g) = ff2(t)*^'(g) 



(7) 
(8) 



.T.(l)/ 



The spatial part * (q) is the initial displacement field (this is 
Zel'dovich (1970) solution), while the second order obeys 



v.*^^^ = E 



8<J 



5*(') d¥p 5*^'' 5*$" 



dqi dqj 



dqj dqj 



(9) 



2.2. Redshift distortion in the Lagrangian framework 

Let us now consider the case of spherical coordinates, when 
radial distances to the observer are estimated by means of red- 
shift measurements. In that case, peculiar velocities distort the 
clustering pattern in redshift space, and the large scale infall 
increases the apparent density contrast. We will now use the 
Lagrangian formalism to derive an analytic expression, relat- 
ing real space and redshift measurements, including weakly 
non-linear corrections. 

The redshift space comoving position z of a particle located 
in r{q) = a x(q) is 



da\ 
'dt ) 



dr r 

dt r 



(14) 



Its modulus is then given at second order by 

z = qr+s[\ + fdt)]9im'-r\q) 
+ S^[l + f2(t)]92m?(q)+0(s'), 



(15) 



where qr stands for q ■ r/r. In the above expression, we have 
explicitly used the separability of *'^'' and ^^ \ and defined 
the logarithmic derivatives of the growth rates 



/i = (a/ 91) dgi/da, and /2 = (a/ §2) dg2/da. 



(16) 
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Limited expansions of /i and /2 near il = 1 yield, in the A = 
case, /i X il'*/^ and /2 x 20^/^. On the other hand, a better 
analytical approximation for il between 0. 1 and 1 is given by 



/i ^ fi^/^ and /2 f« 2 0^*/^ when A = 0. 



(17) 



The approximation for /i was proposed by Peebles (1976), 
while BJCP proposed the one for /2. For a flat universe with 

A j^ 0, BCHJ suggest 

/l«fi5/^ and /2«2fi'>/", whenfi + ^ = l. (18) 

This /i term is very close to the fits previously proposed by 
Lahav et al. (1991) and Martel (1991). 

With the coordinate transformation (15), we could directly 
compute the density contrast in redshift space according to the 
mass conservation requirement, that reads. 



6z(z(q)) = 



dz 



dq 



1. 



(19) 



Instead, we proceed in two steps: we map the Lagrangian q 
coordinates onto the Eulerian ones x, and then proceed to map 
the latter onto our z coordinates. We then see explicitly the terms 
brought by the real-to-redshift space mapping. Additionally, it 
makes it easier to handle the fact that the vector q is in general 
not parallel to t- , or in other words qr "f q- The mapping from real 
space to redshift space is then simplified by using the Eulerian 
coordinate x = r/a as an intermediate step in the calculation, 
because z oc r. After some algebra (see Appendix A.l), the 
density contrast in redshift space can be written as 



SAz) = s{Si + Al) + sHSi + A2) + 0{s^), 



with 



(1) 



6i = -V,.^ 



62 = V^.* +* .V^(V^.* ) 



and 



1 + ^ 1 V,.*^ 

92, 






(20) 

(21) 

(22) 
(23) 

(24) 



A2 = /i 



z'^oz z^oz 



d 



'-.m. 



oz 



+ /? 



' iz^¥P)U^P .3^^ 



z^dz 



dz 



+#(1) A I ^(^2^(1)) 



dz I z'^dz 



d 



z^oz 



(25) 



Note the introduction of the new displacement fields $ and 

" (2) 

* defined by 

*^"(t, z) = *("(t, q = z), ^^^\t, z) = ^^^\t, q = z). (26) 



By definition, they satisfy Eqs. (9), (10) and (11), but with all 
derivatives taken and calculations made in z space instead of q 
space. In real space, i.e., when one is interested in the calcula- 
tion of b{x), the terms Ai and A2 vanish, which is equivalent 

to setting f\ = fi = 0, and the displacements ^ (t, z) and 

" (2) 

* (t, z) must be replaced by the Eulerian displacement fields 
*^"(t, x) = *<"(t, q = x) and *^^'(t, x) = *<^'(t, q = x). 

The expressions Ai and A2 contain terms of the form 
F{z)/z or F{z)/z^, where F(z) is a quantity without preferred 
symmetry nor direction, because we assume the universe is spa- 
tially homogeneous and isotropic. When z is large compared to 
the scale £ considered (which can be for instance a smoothing 
scale, a separation between two objects, or simply a Fourier 
space wavelength), these terms can be neglected compared to 
those in F{z). From a statistical point of view, this requires the 
volume sampled, V, to have a size much larger than £. In that 
case, when statistical quantities are measured in that volume 
V, even if the small z contributions are taken into account, the 
terms in {F(z)/z) and (^F{z)/ z^'^ have a negligible weight. We 
now assume this "large volume" limit. The quantities Ai and 
A2 then simplify dramatically and read 

Ai = -/i^*^", (27) 



A2 



'dz " ' 

/i7^{V,.(* 

oz 

oz 



oz 



(1)#('\t , .2 



* 



(1) 



dz 



(28) 



Our assumption is almost equivalent to the "infinitely remote 
observer" approximation used by BCHJ. But there is a sig- 
nificant difference: instead of using a Cartesian grid deformed 
along the line-of-sight, we keep the radial feature of the redshift 
distortion, i.e., we remain in spherical coordinates. 

3. Redshift distortion of the bispectrum 

In this section, we use a statistical approach and work in 
Fourier space. We first compute the bispectrum5(fci, fc2, fcs) = 
{^ki^k2^k2)k +k,+k =0 i"^ '"^'l^hift space, in the weakly non-linear 
regime. Then we quantify the effect of the redshift distortion 
by using the reduced amplitude Q of the bispectrum which was 
introduced by Fry & Seldner in 1982 (hereafter FS). 

The initial fluctuations are assumed to be Gaussian- 
distributed. In Fourier space, this is equivalent to 



(29) 
(30) 



6i(k)6i(l)) = 6uirAk + l)Pi{k), 



6i(k)6i(l)6i(m) 



0, 



6i(k)6i(l)6i(m)6i(n) ) = SuirAk + 1) SDirAm + n) 

Pi(k)Pi(m) + perm. (31) 
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and so on. Here "perm." stands for the terms obtained by pair 
permutations of k, I, m,n (3 terms in all). The function Si(k) is 
the Fourier transform of the initial density contrast Si{x), and 
Pi{k) =< \Si(k)\^ > is its power spectrum. 

It is convenient to evaluate the quantities Si, S2, Ai and A2 
defined in Sect. 2.2 by using Fourier analysis. The correspond- 
ing integrals are given in Appendix A.2 [Eqs. (A16), (A18), 
(A 19), (A20)]. After taking their ensemble averages, we get at 
lowest order in s. 



Again, we can set /i = /i = to have the value of the bispectrum 
in real space and get the well known expression (obtained by 
Fry in 1984 for fi = 1) 



B{kuk2,k3) = P(ki)P{k2) 



P,(k) = hP(k), 

with 

2 1 . 



(32) 



(33) 



This is nothing else but the result of Kaiser [1977, Eq. (1)], 
when averaged over angles. The bispectrum in redshift space is 
given by 

B^{k\, fc2 ,kT,) = P{k\)P{k2)[ao + 02 cos, 6x2 + 04 cos, 6\2 

( n 3 /) \ /^^l ^2 

+ (fl; COS U\2 + as COS U\2) -, H -— 

\K2 Ki 



+ (1 - COS^6'i2 

+ eye. 



2 /x/ {k\ + k'^ +X2k\k2COs6\2 
\ k\ + k\ + 2k\k2Cos6\2 



(34) 



where cos 6ij = ki.kj /(kikj). Here "eye." stands for the terms 
obtained by cyclic permutation of the triple-valued indices. The 
factors ai, i = 0, . . . ,4, xi and X2 depend on the values of il 
and A. We list them here: 



ao = (l-'^) + |(2-«;)/i + ^(ll-«;)/f+^/? + ^/f 



^3^ 15^ 



«;/2(^ + fi/i + iifh 



35' 



ai 



, 4„ 14 „2 12 „3 15 „4 

1 + -fi + — ft + — n + ft, 

3 15 ' 35 ' 315 ' 



02 = (1 + «:)+ =(2 + «:)/! + — (23 -«:)/? +^/,^ + ^/f 



I' 

1 _6_ 

'3 ■ I5' 



15^ 



30 

35' 



48 
315' 



+ «;/2(- + 77/1 - i^fh 



3_ 
35- 



"^ 15 ' "^ 35 ' "^ 315 " 



Xl = 

X2 = 

and 



15 



Kf\f2 



-^«/?/2, 



K = gi/g\- 



(35) 



(36) 



1.^1(1 
9\ 



2H-COS 6112 

COS^ 6*12) 



fcl_ ^2 

k2 ki 



H-cyc. 



(37) 



The values of ai/b], Xj/b^ are listed in Table 1 for 



(fi,A/3i?2) = (i^0),(0.1,0)and(0.1,0.9).Exceptwhenfi= 1, 
we have (ai/b^)^ ~ (ai/&i)reai- In other words, when il is suf- 
ficiently small compared to unity, and A is arbitrary, we obtain 

BAkuk2,k3)-bJB(kuk2,k3). (38) 



When il approaches unity, the differences between real space 
and redshift space are mainly introduced by the factors ao/bj, 
a2/&i and a^/b^. 



Table 1. Factors 61, a, and Xj involved in the computation of the 
power-spectrum and the bispectrum for various values of Q 





1 


■edshift space 




real space 




Fact. 


Q = l 


Q = .l 


= .!" 


Q = 1 


Q = .l 


= .!" 


61 


1.867 


1.180 


1.201 


1 


1 


1 


ao/bi 


1.335 


1.459 


1.429 


1.429 


1.461 


1.436 


ai/61 


1.050 


1.005 


1.006 


1 


1 


1 


a2/bj 


1.043 


.5770 


.6145 


.5714 


.5389 


.5644 


a,lh] 


.1603 


.0149 


.0180 











aAlh] 


.0411 


.0034 


.0044 











-x,lh\ 


.0328 


.0060 


.0064 











xilh\ 


.0281 


.0013 


.0015 












" For a flat Universe with A/3-ff ^ = .9 

To illustrate our results, it is convenient to use the reduced 
amplitude Q of the bispectrum. It is defined by (FS) 



Q(kuk2,k3) 



B(ki,k2,k3) 



P{k,)P{k2) + P(k2)P(k3) + P(k3)P(ki) 



.(39) 



The parameter Q, defined above, must satisfy Eq. (34) and the 
constraint ki + k2 + k3 = 0. Therefore, Q can be described as 
a function of ki , k2 and cos 6 = cos 6112. It also depends on the 
shape of the initial power spectrum and was first measured in the 
galaxy distribution by FS. Comparisons of theoretical predic- 
tions with #-body experiments were done by Fry et al. (1993). 
In the weakly non-linear regime, they found a good agreement 
between the numerical measurements and the predictions of 
second order perturbation theory. 

Figure 1 gives Q as a function of 6 for ki/k2 = 2 and 
various values of (fi, A) [(1,0), (0.1, 0) and (0.1, 0.9)]. We have 
considered 5 scale-invariant power- spectra, with n = —3, -2, 
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Fig. 1. The reduced ampUtude Q of the bispectrum in real space versus 
the wavevectors angle 6 in the case k\/k2 = 2, for Q = 1 (solid line), 
Q = 0.1 and A = O(short dashes) andQ = 0.1 andQ+A/3-ff^ = 1 (long 
dashes). Scale-Invariant power spectra are considered with ra = — 3 , 
-2, -1, 0, +1 (from top to bottom). 



-1, and +1. In that case, Q = Q{ki/k2,cos6i2,n,i2). As 
expected, the dependence of Q on il and A is very weak. The 
first panel of Fig. 2 is the redshift space analog of Fig. 1 . We have 
kept the solid curves, which give Q in real space for il = 1 . As 
we already saw above, there is a significant difference between 
the real and redshift space values of Q only when il ^ 1 . This 
difference can be as large as ~ 20%, but the general features 
of Q versus 6 and ki/k2 are not changed, as confirmed by the 
second and third panels of Fig. 2 which shows Q versus 6 for 
ki/ki = 1 and for ki/k2 = 10 in real and redshift space when 
fi= 1. 

Finally, when ki/k2 = 2, the ratio Qz/Q depends only 
weakly on the initial power spectrum, as can be seen in Fig. 3. 
We find that the il-dependence of this ratio for il ^ 1 is given 
by 



Qz(ki,k2, kj) = Q(ki,k2, k3)[l + A{ki, ^2, ^3)], 



(40) 



where A{ki , k2, kj,) is a smoothly varying function such that 



\A{ki,k2,k3)\^Q.2n. 



(41) 



So far, all of the quantities derived here were statistical measures 
of the spatial distribution of matter. However, if the galaxy 
distribution is biased with respect to that of matter, the reduced 
moments (like the skewness), estimated from counts of galaxies 
do not have to coincide with those for the matter. Following Fry 
(1994), we will now rederive the effect of a local biasing on 
the Q{k\/k2, 0, i2, n) parameter. The local biasing prescription 
assumes that the number-density field of galaxies, ng(x) is a 
smooth local function / of the underlying density field of the 
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Fig. 2. The reduced amplitude Q of the bispectrum as a function of 
the wavevectors angle 6, in real and redshift space for several values 
ofki/k2. As in Fig. 1, we consider scale-invariantpower-spectra with 
n = —3 (top curve), -2, -1, 0, +1 (bottom curve). In all panels, the full 
line and the dotted line refer respectively to real and redshift space for 
Q = 1. In the top panel, the short dashes correspond to Q = 0. 1, A = 
and the long dashes to Q = 0. 1, fi + A/3H^ = 1, in redshift space. 
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matter, p(x), that is Ugix) 
expansion 



6g = b6 + T,'^^2- 



hmS" 



f [p(x)]- This allows a Taylor 



(42) 



where b is the usual "linear bias", while higher order terms allow 
for nonlinear biasing. They are also necessary to ensure that 
each term in the perturbative expansion for Sg has a zero mean 
value (cf. Fry & Gaztanaga 1993; Juszkiewicz et al. 1993b). 
To derive the relationship between Q and its estimate based 
on galaxy counts, Qg, it is enough to substitute the expansion 
for Sg(S) into the perturbative series for S, then rederive the 
bispectrum and the Q parameter. The result of all this is simply 
given by (Fry 1994) 



Qs 



Q &2 



(43) 



Thus both the linear bias b and the non-linear bias &2 can be 
determined by comparing the measured shape of Qg{k\/k2, 0) 
with the theoretical predictions for the underlying matter field 
(assuming that indeed large scale structures arose from the grav- 
itational development of Gaussian initial conditions, which can 
be tested otherwise). This may solve a well known problem 
which always arises in conventional dynamical determinations 
of the mean density of the universe. Such studies tend to produce 
only estimates of il^^/fe, where the linear bias parameter, b, is 
unknown (see, e.g. Strauss et al. 1992). As a result, a biased 
il = 1 model is dynamically indistinguishable from an open 
unbiased one. Therefore, in conjunction with our calculations, 
which allow to properly take into account the effects of redshift 
distortion on Q, Fry's result may for the first time provide a way 
to use redshift surveys to break the degeneracy, and decouple 



the estimates of linear and non-linear bias from the estimates 
of il and A. 

4. Skewness of the smoothed density field in redshift space 

We now leave Fourier space and turn to the skewness factor 
53 = [S^) I ip^') of the smoothed density field. The smooth- 
ing windows considered here are the most widely used ones, 
i.e. the Gaussian window and the top-hat window. Firstly, we 
study analytically the redshift space distortion in the weakly 
non-linear regime. Then, we compare analytical results with 
careful measurements in #-body experiments and estimate a 
dynamical range of validity of these results. 

4.1. Analytical predictions in the weakly non-linear regime 

In the weakly non-linear regime, one can use second order per- 
turbation theory to compute 53 (JB, JBC and Appendix A.2). 
By using Fourier integrals, the low orders moments of the den- 
sity field ((5^) and 1^8^') can be written in terms of, respectively, 
the power spectrum and the bispectrum. This gives 

a\l)={8^)= [J^P(k)Wl(k), 



(27r) 



(44) 



6') = 



dkjdky 



■ B{ki,k2,-ki - fc2) 



(27r)6 

W^(ki)W^(k2)W^(\ki + k2\), 



(45) 



where W(,(k) stands for the Fourier transform of the smoothing 
filter of characteristic size £. These expressions are of course 
also valid in redshift space. But then the quantities P(k) and 
B{ki , k2, kj,) have to be replaced by their redshift space coun- 
terparts Pz{k) and Bz(ki,k2,kj). Appendix B gives a full 
length account of the computation of the redshift space skew- 
ness (and hence of the real space if one takes /i = /i = 0) for 
scale-invariant initial power- spectra, with or without smooth- 
ing. Here we just give the main results. We first consider the case 
without smoothing. In some cases, the integrals in Eqs. (44), 
(45) may diverge. In order to insure convergence, one then 
needs to introduce a large-A; cutoff in the power spectrum. Nev- 
ertheless, the ratio 53 = [S^) / [S^) may still be well defined, 
provided it remains finite when the cutoff tends to infinity. In 
redshift space, it is given by 



5- 



3,z 



3 / 02 04 ~ 

77 [ao+ — + — +X 



(46) 



where X depends on the shape of the power spectrum, and is 
hard to compute analytically. Nevertheless, it can usefully be 
bounded as follows 



29 23 - 19 7 

X] X2 < X < X] H X2. 

20 120 20 120 



(47) 



By taking /i = /2 = in Eq. (46), (35), and (36) we recover 
the familiar value of ^3 in real space, that is 

,92 



53 = 4 - 2- 



9t 



(48) 
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which is independent of the initial power spectrum and nearly 
independent of il, as found by BJCP and JBC. The following 
values of 53 ^ , 



Table 3. Skewness in redshift and real space for different power spectra 
P(k) = k" with top-hat smoothing 



fi= 1, A = 0, 

fi= .1, A = 0, 

fi = .l, A/3H^=.9, 



S3,z =5.01 ±0.02, 53=4.86, 

53,^ = 4.95 ± 0.01 , 53 = 4.92, (49) 

Si, =4.89 ±0.01, 53=4.87, 



show that the projection in redshift space does not significantly 
change the value of the skewness factor of the unsmoothed 
density field. This confirms the results of BJCP and BCHJ. 
Indeed, by using their "infinitely remote observer" approxima- 
tion, BJCP and BCHJ had found 83^, = 5.07 ± 0.07 for fi = 1 
and 53,^ = 4.96 ± 0.03, for fi = 0.1 (A = 0). These results are 
in a good agreement with those obtained here, as expected. 

Numerical values of S3 in real and redshift space for a 
Gaussian and a top-hat smoothing are listed in Tables 2 and 3. 
Figure 4 summarizes the behavior of skewness when the spec- 
tral index n is varied, for background model with A = 0, and 
il = 0.1oril = l.JB and JBC discussed rather extensively the n 
dependence of S3, and JBC had found that S3 is nearly indepen- 
dent of il. Here we note that this remains true for spatially-flat 
universes with non-vanishing cosmological constant. 

The mapping from real to redshift space brings only very 
small changes, except for values of il close to unity. In the case 
il = 1, the difference between 53 _^ and S3 increases with n, but 
remains small, at most of order 15%. This behavior is maybe 
slightly different from what one could have naively expected 
from the analysis of the bispectrum (e.g. Fig. 3). In that case, 
indeed, the effect of a redshift space projection appeared to be 
of the same order for any value of n. For a top-hat filter, the 
linear relation between S3 and n in real space maps into another 
linear relation between Si^z and n. For example, when il = 1, 
S3^z is well fitted by 



S3,, 



35.2 



1.15(rn-3). 



(50) 



This is to be compared with the real space values given by 
Eq. (3). 



Table 2. Skewness in redshift and real space for different power spectra 
P(k) = k" with Gaussian smoothing 









Skewness 
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redshift 


space 




in real space 
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Q = l 


Q = .l 


= .!" 


Q = l 


Q = .l 


Q = .r 


-3 


5.022 


4.947 


4.894 


4.857 


4.922 


4.871 


-2 


4.043 


4.1 


4.042 


4.022 


4.089 


4.036 


-1 


3.364 
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3.475 
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3.541 


3.484 





2.914 


3.217 


3.139 


3.144 


3.228 


3.162 
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2.652 
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3.012 


3.029 


3.132 


3.051 



For a flat Universe with A/3H^ 









Skewness 








in 


redshift 


space 




in real space 


n 


Q = 1 


Q = .l 


= .!" 


Q= 1 


Q = .l 


= .!" 


-3 


5.022 


4.947 


4.894 


4.857 


4.922 


4.871 


-5/2 


4.443 


4.439 


4.384 


4.357 


4.422 


4.371 


-2 


3.868 


3.932 


3.876 


3.857 


3.922 


3.871 


-3/2 


3.294 


3.425 


3.367 


3.357 


3.422 


3.371 


-1 


2.72 


2.917 


2.858 


2.857 


2.922 


2.871 


-1/2 


2.144 


2.41 


2.349 


2.357 


2.422 


2.371 





1.557 


1.902 


1.839 


1.857 


1.922 


1.871 



" For a flat Universe with A/3H^ 



4.2. Measurements of S3 in il = 1, scale-invariant, N-body 
simulations 

Numerous previous comparisons with #-body simulations 
showed a very good agreement with second order perturbation 
theory predictions up to a (density contrast) variance of order 
unity (see, e.g., JBC and Juszkiewicz et al. 1993b). However, all 
these measurements were made in real space. Here we extend 
these tests to the predictions concerning the redshift distortion 
in order to determine their range of validity. In particular, it is 
important to find out, when strongly non-linear effects, like the 
"fingers of god", become large enough to invalidate our pertur- 
bative results. Section 4.2.1 describes the simulations done, the 
measurements method, and our error estimates, while section 
4.2.2 discusses the results. 

4.2.1. Measurement: method and uncertainties 

Using the PM code of Moutarde et al. (1991), we performed 
three #-body simulations with il = 1 , 64-' matter particles in 
a 128^ grid, and scale-invariant initial conditions with P(k) = 
Ak",n = 0, -1, -2. Following Efstathiou et al. (1988), we 
picked a normalization of the initial power spectrum so that it 
matched the white-noise level at the Nyquist frequency of the 
particle grid, except for n = —2. In that case, the normalization 
constant A was chosen 4 times smaller. The density field was 
estimated by convolving the particle distribution with either a 
Gaussian filter of half- width £, or a spherical window of radius 
£ (all lengths are expressed in units of the size of the simulation 
box Lbox = !)■ 

In principle, the only relevant scale in the simulations is the 
correlation length £0, for which the two-body correlation func- 
tion ^2 is unity. Thus measurements made at different times 
of a simulation should give statistically equivalent results pro- 
vided there are made at the same fixed fraction of that (time- 
dependent) scale £(). As far as S3 is concerned, it should thus 
be constant at constant variance (for a given value of n). If the 
expected scaling behavior is only approximate, it may come 
from one or several of the following numerical limitations: 
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Fig. 4. Slcewness factor in real and redshift space, for Gaussian and top-hat smoothing as a function of the power spectrum index n. 



Initial conditions: they were created by slightly displac- 
ing matter particles from a regular pattern by using the 
Zel'dovich approximation. This is not totally innocuous. 
Firstly, the constraint arising from this regular pattern may 
contaminate the latter measurements (in particular in under- 
dense regions, see e.g. Bouchet et al. 1991, Bouchet & Hern- 
quist 1992, Colombi et al. 1994). Secondly, the Zel'dovich 
prescription is only valid at first (linear) order, and thus 
does not set up properly second order quantities such as the 
skewness (JBC). One thus has to wait long enough for tran- 
sients to become negligible. Typically, the correlation length 
needs to be larger than the mean interparticle distance. This 
means that we should have 



^ 1/64 ~ 0.0156. 



(51) 



Finite spatial resolution: Bouchet et al. (1985) showed that, 
in the linear regime, the effect of the grid could affect the 
linear growth of fluctuations on scales substantially larger 
than the grid scale Ax (say 8Aa; for a typical PM code 
without a staggered mesh, if better than 10% accuracy per 
mode is required). On the other hand, the transfer of power 
in the gravitational clustering problem appears to proceed 
from the large to the small scale. This means that "gravity is 
forgiving", and reasonably accurate results in the strongly 
non-linear regime can be obtained even on scales compa- 
rable to the cell size (e.g. Weinberg et al. 2001, but PM ve- 
locities may systematically be underestimated by as much 
as 20%). Thus provided that our measurements are such 
that the correlation length £() remains large as compared 
to the grid scale (but small as compared to the box size), 
our skewness measurements should be accurate. Indeed, 
this insures that measurements on scales comparable to the 
grid scale (logjgf ^ — 2.1) are in the strongly non-linear 
regime, and the weakly non-linear measurements are on 
much larger scales, where a PM code is reasonnably ac- 
curate. Furthermore, we compared our results to those ob- 
tained from high resolution simulations (made with the tree 
code of Bouchet & Hernquist, 1988) with identical initial 



conditions (Colombi et al. 1994). The differences are rather 
small, of order at most ten percent; the measured 53 in the 
high resolution simulations are only slightly smaller than 
the measured 53 in our PM simulations. 

3. Finite volume effects: since we used periodic boundary con- 
ditions, the power coming from scales larger than the box 
size is missing. Also, at scales smaller but close to the box 
size, only a few independent modes of the power spectrum 
are sampled. This increases the uncertainties on the mea- 
surements as one approaches the box size. To minimize 
these effects, one usually imposes logjQf ^ — 1. But this 
constraint is not necessarily sufficient. Indeed, the large p 
tail of the PDF is determined by just a few large clusters and 
that tail is thus subject to fluctuations due to small number 
statistics, till it reaches an artificial cutoff at pmax (Colombi 
et al. 1994, hereafter CBS). If the volume is large enough, 
these modifications of the true tail happen for such small 
values of the PDF that is has little effect (however, higher 
order statistics are more sensitive to this scale-dependent 
effect than those of lower order). If the correlation length, 
which gives the typical size of a cluster, is small compared 
to the box size, then this effect should also be small and can 
be corrected for (or at least the corresponding error-bar can 
be evaluated, see CBS and Appendix C). 

4. Discreteness effects: starting from a discrete particle repre- 
sentation, we want to measure statistical properties of the 
underlying continuous density field. This shot noise can eas- 
ily be corrected for a top-hat filter and is naturally erased 
by a Gaussian filter. It can be neglected at scales large com- 
pared to the typical distance between two particles in a 
cluster (e.g., Balian & Schaeffer 1989). 

The considerations above lead us to measure a^ and 53 in 
the following intervals 



-2.1<logioi'<-1.3, Alogioi' = 0.2, 
for a Gaussian filter and 



2.0< logio£< -1.2, A log 



10' 



0.2, 



(52) 



(53) 
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Fig. 5. Skewness factor in real space (filled circles) and redshift space (open circles) measured in TV -body simulations. The dashes show the 
theoretically expected values in real space and the dots-dashes to the redshift space one, as given in Tables 2 and 3. The left panels correspond 
to a Gaussian smoothing and the right panels to a top-hat smoothing.The initial power-spectra considered here are scale invariant, with ra = — 2 
(top), -1 (middle), and (bottom). 
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for a top-hat filter. The quantity A logjgi' corresponds to the 
logarithmic scale step we chose. We used several snapshots of 
each simulation, corresponding to the various expansion factors 
a (with initially a = 1) given in Table 4. For each snapshot, the 
constraint (51) is approximately fulfilled and £q remains rea- 
sonably small (as compared to Lbox, Table 4). With such values 
of £() we expect finite volume effects to be small and we do not 
correct for them. However, in order to get an idea of the uncer- 
tainties related to the possible artifacts listed above, we tried 
to compute as realistic error-bars as possible. In particular, we 
used the procedure proposed by CBS to estimate finite volume 
effect errors (see Appendix C, where we discuss in detail the 
error-bars calculations). 



Table 4. The expansion factor a and the measured correlation length £q 
corresponding to the various snapshots chosen to measure the skewness 
in our TV -body simulations 



n 


a 


£o 




2.5 


0.027 


-2 


3.2 


0.040 




4.0 


0.054 




2.5 


0.015 


-1 


4.0 


0.022 




6.4 


0.036 




4.1 


0.011 





8.1 


0.018 




16 


0.028 



4.2.2. Measurement results 

Figure 5 displays the measured values of 53 in real (filled cir- 
cles) and redshift space (open symbols) as functions of the 
logarithm of the variance. Left panels correspond to Gaussian 
smoothing and right panels to top-hat smoothing. The index n 
of the initial power spectrum increases from top to bottom. The 
dashed and dotted-dashed lines represent the analytical predic- 
tions of second order perturbation theory in real and redshift 
space respectively. They should superpose to the filled circles 
and the open circles in the weakly non-linear regime, i.e., when 
<t2<1. 

In real space, there is a good agreement between theory 
and measurement in the regime a ^ I, which confirms the 
results of previous authors. In redshift space, the agreement 
between the second order predictions and the measurement is 
globally less good (except for the case n = — 1 with top-hat 
smoothing), unless a^ is much smaller than unity, typically 
(T^ ^ 0.1. Indeed, the measured quantity Sj^z is seen to be 
considerably smaller than 53, particularly in the strongly non- 
linear regime. The Ss^z vs. (T^ curve tends also to be flatter than 
its real space image, 53(17^), as already noticed by Lahav et al. 
(1993); the overall deviation from a scale-invariant behavior 
[53(17^) = constant] is also much less significant in redshift 
space than in real space. This deviation, which seems in real 
space to increase with decreasing n (in agreement with the 



results of Lucchin et al. 1994) is particularly large for n = —2. 
In the latter case, the measured value of S3 is amplified by more 
than a factor two as compared to the value seen in the weakly 
non-linear regime (left lower corners of top panels in Figure 5). 
Note however that the case n = —2 is expected to be strongly 
contaminated by finite volume effects (as already discussed by 
CBS and Lucchin et al. 1994). This is why the error-bars are so 
large in that case. Furthermore, the measurements of S3 and S3^z 
are likely to be underestimates of the true value (See Appendix 
C and CBS). 

The fact that S3^z is much smaller and flatter than S3 in 
the non linear regime is certainly related to the "finger of god" 
effect described in the introduction. This is of course a highly 
non-linear effect. It has been analytically studied for the three- 
body correlation function by Matsubara (1994). His results, 
when generalized to averaged correlations should match our 
measurements. Furthermore, it is likely that this finger of god 
effect is also responsible for the slow convergence toward the 
perturbative predictions when the variance is decreased. Indeed, 
it smears dense structures along the line-of-sight up to fairly 
large scales. This behavior has already been observed in stud- 
ies of redshift distortions of the power spectrum (e.g. Cole et 
al. 1994). 

5. Conclusions 

In this paper, we have studied redshift distortions of clustering 
in the weakly non-linear regime of the gravitational instability. 
Second-order Lagrangian perturbation theory was applied to 
spatially flat cosmological models with il + A/3H^ = 1 and to 
open models with A = 0. We have analytically estimated the 
bispectrum and the redshift space value S3^z, of the skewness 
factor 53 = (S^) / (S^) for both the unsmoothed and smoothed 
density field. We considered the case of Gaussian and top-hat 
smoothing for scale-invariant initial power-spectra (|i5fcp) oc 
k" . We also compared our analytical results to measurements 
in #-body experiments (for il = 1). The main results are the 
following: 

1 . The skewness factor and the Q-amplitude of the bispectrum 
are seen to be relatively insensitive to projection in redshift 
space in the weakly non-linear regime. In other words, 
the third moment and the square of the second moment are 
globally changed by nearly the same factor [~ fe^Eq. (33)], 
both in real space and redshift space. The most important 
differences are found for il = 1 , where | S'3_^ — 53 1 /5'3 ^ 0.2, 
which is small in view of the present uncertainties in the 
observationnal determinations of S3. 

2. The skewness, which was already known to be almost inde- 
pendent of the value of the density parameter il (for A = 0), 
depends only very weakly on the value of A in a flat uni- 
verse (il + A/3H^ = 1). The same result is valid for the 
bispectrum. 

3. Analytical predictions for S3 are in good agreement with 
measurements in #-body simulations but the finger of god 
effect is felt till quite large scales . Indeed, our measurements 
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suggest that the #-body values of Sj^z do converge to the 
analytical predictions only in the regime a^ ^ 0.1 (with 
a^ = (S^))- In real space, the agreement between theory 
and measurements was quite good for a almost as large as 
unity. 

Peebles (1980) had shown that 53 for the unsmoothed field 
has a specific value, 34/7, if large scale structures arose from 
the gravitational development of Gaussian initial conditions 
in an il = 1, A = 0, model. Bouchet et al. (1992) showed 
that the effect of il and of the real space-redshift space map- 
ping are quite weak. Juszkiewicz and Bouchet (1991) and 
Juszkiewicz et al. (1993a) computed the power spectrum de- 
pendence introduced by the smoothing of the evolved density 
field. Juszkiewicz et al. (1993b), and Fry and Gaztafiaga (1993) 
analyzed the effect of a non-linear biasing. This paper com- 
pletes the series by an analysis of the redshift space value of 53 
for a smoothed field in a large class of cosmological models. At 
all steps the results were carefully checked with simulations. 
Comparisons with redshift space measurements like those in 
the IRAS 1.2 Jy catalog (Bouchet et al. 1991, 1993) can then 
be used to test the hypothesis that the presently observed large 
scale structure was generated by gravitational instability, acting 
on initially Gaussian density fluctuations. 

The real space expression for the bispectrum was first com- 
puted by Fry (1984), while the region of validity of this per- 
turbative result was tested in numerical simulations by Fry et 
al. (1993). Most recently. Fry (1994) pointed out that measure- 
ments of the bispectrum can be used to measure the value of 
the non-linear bias. Here also we complete the series by ana- 
lyzing how il and A affect the bispectrum and the Q parameter, 
measured in redshift space. As a result, it may for the first time 
become possible to use redshift surveys to estimate the linear 
and non-linear bias factors, breaking the coupling between the 
bias and cosmological parameters - il and A. 
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APPENDIX 

A. From real space to redshift space 

Let us consider a spherical coordinate system centered on the 
observer. The mapping from real space to redshift space can 
then be written as 



ax = (r, 9, Lp) 



■\z,9,^), 



(Al) 



where z is the "comoving" redshift given by z = Vr/{da/dt) 
and Vr is the proper velocity projected along the line-of-sight. 
In the Lagrangian approach, the comoving coordinate of a fluid 
element (a "particle") is 



x(t) = q + £*^"(g, t) + £^*^ '(g, t) + C(£ ) 



2.7,(2)/- 



(A2) 



We assume here that the fastest growing modes, gi{t), i = 1,2, 
are dominating. We can thus write ^^^'{q, t) = gi(t)^ (q, t). 
The comoving redshift coordinate of the element of matter is 
then 



z(q) = q^+e(l + f,)^'^\q,t) 



(A3) 



with fi = {a/ gi){dgi/da). We also denote by a subscript the 
radial component of a vector, i.e. u^ stands for the radial com- 
ponent of the vector u,Ux = u ■ x/x = u ■ r/r. 

By introducing the new vectors ^ (x,t) = ^^^\q = x,t) 
and ^^ \x,t) = ^''^\q = x,t) which obey equations (9), (10), 
but in X space instead of q space, we can express z in Eulerian 
coordinates as 



z{x) = x + efi'^^^Xx) 



+ s^f2^'i\x)-fl^'■'\vM^)+0(s') 



(A4) 



A.l. density contrast 



We now write the density contrast in redshift space, as a function 
of z. In redshift space, the density contrast S^ = (pz — p)/p 
around a particle initially in q is given by the requirement of 
mass conservation: 



6z(q) = 



dz 



dq 



1. 



(A5) 



The direct use of Eq. (A5) is however a little awkward, because 
the vectors z and q are not necessarily collinear. We found 
simpler and more illuminating to first compute the density con- 
trast in real space as a function of the Eulerian coordinate x, 
using mass conservation requirement (5) and then pass to red- 
shift space using the transformation (A4). The development at 
second order in e of Eq. (5) is 

8Aq) = -£V„.*(" 



+ s 



(V< 



^^mf + 



1 + ^ 

91 



V,.*<2> 



(A6) 



v(«)/ 



Using the Eulerian displacements * {x^f) defined above and 
the following property 



V,.*(" = V,.*^" - e {v,(V,.*^")} * 



^" + 0(e^), 



(A7) 



we can easily get b^ as a function of x at second order in e 
(e.g. Bouchet et al. 1994). Then, using the fact that z ex a;, we 
simply write the mapping from real space to redshift space as 



Px 



dx 



dz 



(A8) 
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By using Eq. (A4), we are now able to write the density contrast Since 
in redshift space as a function of x : 



-.Wt 



pAx) = l-£(V.*^" + /i^(a;2^(,i)) 



+ £ 



+/l 
-fl 



x^dx 



92 



6i = -V,.*^"(z) = -gi(t)J2^i(ky 

k 

we get 



ik.z 



x^dx ^ x^dx 

r) r) \I/(1)2 



^k 

and 



k 



6i(k), 



x^dx 
d 



dx 



x^ 



,2 ik.z 



x^dr 



(x^^ff^) 



(A9) 



k 

Hz,t) = gj(t)Y,'^i(k)6i(l)W(\k + l\) 



(A16) 



(A17) 



(A18) 



The last step is to express S^ as a function of z. We introduce 
the new vectors * (z,t) = ^'^ (x = z,t) and # (z,t) = 

- (2) 

* (x = z,t)- These displacements obey the equations (9) 
and (10), but with the calculations done in z space. Using the 
following properties, 



k,l 



1 + COS ( 



gf 2 



i(k+l).z 



(A19) 



V..*<" = V..#^'' 



2,Tr(Ih 



d 



.(IX 



x'^dx 



(a;^*^'0 



z^dz 



sfi¥P—(v,.¥'), 

oz 



{zWP) 



(AlO) 



A2(z) = g\{t)Y,'6Ak)lx{l)W{\k + l\){kiik + lpi) 

k,l 
/ly(l+COs6'y)+/fy//f 

gz l -COS^ 6* (fc/jfc +///;) ^ gi(fc+0^^ 



9l 



' dz [z^dz^ ' ' 



(All) 



with 



Pk 



k.z 

kz 



^ii 



Iz 



COS! 



(fc + rp- 
_ k.l 



(A20) 



(^?^2) 



In this paper, we assume that the linear density contrast, 
we obtain the density contrast in redshift space as a function Si, is Gaussian. Its Fourier coefficients, gA, are thus uncorre- 
of z, as given in Eqs. (20), (21), (23), (24) and (25). In the lated and satisfy Eqs. (29), (30), (31). This implies for instance 
following, we use the "large volume approximation" introduced [Eq. (30)] that the first term of (A13) vanishes. The first non 
in Sect. 2.2, which results in a considerable simplification of the trivial term of (S^) is then 
expression of S^iz) [i.e. Eqs. (27), (28)]. 

The second and third moments of the density contrast are 
computed by ensemble averaging: 

{6Hz))^^ = e^ {8] + 26iAi + A?) (A12) 

{6\z))^ = s'{(Si + Ai)')+3s^S^S2 + S^iA2 

+2(5i Ai(52 + 2(5i Al A2 + AI62 + A^A2) (A13) by using the equations (A13), (A18-A20), (31). In the same 

way, all terms of the second and third moments simplify thanks 
In order to consider the more general case of a density contrast to Eq. (31). The averaging over z of pk and pi factors yields 
field smoothed with a window W(z), 6, (z) in Eq. (20) has to functions of cos 6. Eventually one gets: 
be replaced by 8^{z) = 8^, * W{z) . 



dk dl 



2 + cos ( 



P(k) P(l) W(k) W(l) W(\k + I 
k 



1 + ^ ) (1 - cos^ 6'()4t21) 



A.2. Fourier Transforms 

For a spherically symmetric filter, the smoothed density contrast 
is given by 



■bie^g^i 



6') = 3sV, 



6'^(z) = J26zikWik), 



and the displacement field is 



dk 

(2^ 

dk dl 

(27r)6 



P(k)W\k) 



(A22) 



P(k)P(l)W(k)W(l)W(\k + l\) 



^^'\z,t) = gi(t)J2^k 



) ik.z 



(A14) 



(A15) 



(flO + fl2 COS^ 6+04 cos'* 



+ {ai cos 6 + 03 cos 6) (7 + 7 
+ (1 - cos^ 



k^ 
2 I X] (k^ + fi) + X2 kl cos ( 



k^ + P + 2klcost 



,(A23) 
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where the factors &i, a, and Xj are listed in Sect. 3 [Eqs. (33), 
(35)]. 

Note that by setting /i = /2 = in the expressions above, 
we can recover the real space values of all the quantities just 
computed. JB and JBC had developed techniques to compute 
the skewness of the smoothed density field in real space. These 
techniques will now be described in detail and applied to the 
case of redshift space. The real space results will also be pre- 
sented. 

B. Effect of smoothing 

We now compute the skewness of the density field with either 
no smoothing, a smoothing with a Gaussian filter, or smoothing 
with a top-hat filter. We deal only with power law models, 
P{k) = Ak" . As these models are scale-free, we can set the 
smoothing length to unity without loss of generality. 



flfl H- fl2 M + a4U + [ai u + U} 



+ (1 



X; (a;^ + j/^) + X2 xyu 
x^ + y^ + 2 xyu 




(B4) 



The case n = — 3 is special and must be treated separately from 
the case of the other n. 

B.2.1. Case when n = — 3 

When the power spectrum index n is equal to —3, both ((5^) 
and ((5-') diverge. But we can still give a meaning to 53, if we 
define it as its value in the limit, if any, when n tends to —3. By 
introducing e = n + 3, one gets 

\ / 47r2 4^' 



B.l. Unsmoothed Result 

In the case without any smoothing (W(k) = 1), the integra- 
tion of the equations (A22), (A23) over k and I in spherical 
coordinates yields 



where F is the Euler Gamma function, which behaves as 1/e 
when € goes to zero. For {^^), it is convenient to change vari- 
ables, from the Cartesian {x, y) to the polar ones {p, (p) defined 
hy X = p cos (p,y = p sin (p. The integration over p from to cxd 
then yields 



{s^y 



ao H- — + ^ + ■ 



1/2 



dkk^P{k) 

with u = cos 6. The last term 

(1 -u^)\xi(k^ + f) + X2klu) 



{j^dkk^P(k)) 
f P{1) / duX{k,l,u)\^\) 



X(k,l,u) = 



(k^ + P + 2klu) 



(B2) 



is problematic because k,l, and u cannot be separated. However, 
thanks to the inequalities k^ + fl > 2kl > 0, it is easy to 
bound this term between polynomials of u. Integrating these 
polynomials, we obtain 



29 



23 



-xi 



:X2 < 



' 19 7 

duX(k, I, u) < —xi + T^«2- 



20 120 

Some numerical values are given in Sect. 4. 1 

B.2. Smoothing with a Gaussian filter 



(6') 



3A?_ 
167r4 

cos^~ 

+(1- 



1 .7r/2 

du / d(p 
1 ^0 



r(e) 



(iH-f sin2^)^ 

1 / 2 4 

ip [ao + a2U + a^u 



<^sm 

9^9 {X\ + UX2) ^, -I. X 

M ) -; ; — - — h 2{a\u + uju ) tan ip] 



I + u sin2<^ 
The integrand of Eq. (B5) includes the terms 



(B5) 



(lH--sin2^)" 



1 -elog(lH--sin2^)H-0(0 



which always converges since |M/2sin2<^| < 1, and 



cos 



p sm p 



I +u sin2<^ 



2i-^sin^2^ 

( 

I + u sin2<^ 



H-cos <^sm p. 



In this last equation, the first term on the right hand side is 
bounded as e tends to 0, whereas the second one verifies 



7r/2 



dp cos^ p sin 



.1 ^ r^(e/2) 
^ 2F(e) 



We now consider a density contrast field filtered with a Gaussian 

window of variance unity and normalized to unity. Its Fourier 

transform is then W(k) = exp(—k^/2). The integration over the t 

angular variables of k and I in the equations (A22) and (A23) '''^3('^ = — 3) = -^ 

yields ^ 

, . 00 This yields for instance 

(^2) = £2^ / dx x""-^ cxp-^\ (B3) 

2^^ Jo fi = l, 53 = 5.0217, 



(e.g. Gradshteyn & Ryzhik 3.621.5). Integrals of terms with an 
odd degree in u vanish and one gets for the skewness 



02 04 16 «! 

«°^y^y^i5T 



6') = £^ 



dx 
^0 



dy 



du X 



y exp 



fi = .l, 53 = 4.9469. 
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B.2.2. Case when n is not —3 
The term 



Xix,y,u)= (1 -M ) 



2^ 2 ( xi (x^ + y^) +X2xyu 



x^ + y^ + 2 xyu 



Equation (B6) can be integrated, and one obtains an interval for 
^3 according to the computed bounds on the term X. The three- 
dimensional integral of X is first reduced to a one-dimensional 
integral of an impressively long function which happens to be 
very smooth. This function can then be integrated numerically 
to yield an exact value of the skewness. Numerical values of the 



in the Eq. (B4) makes it impossible to integrate analytically skewness for different power spectra in real and redshift space 
((5^). Fortunately, this term can be bounded. To begin with, are listed m lable 2. 
let us suppose that it is a constant, and let us set m = 2 cos j3. 
Equation (B4) then becomes 



47r4 



27r/3 



dj3 sin j3 



7r/3 



with 

G„(/3) 



(G„(/3) [X +ao + 4a2 cos^/3h- 16 04 cos"^/?) 
+F„(l3) (4fl/ cos l3 + 16 as cos-'/j)) , 



dx / dyx"^^y"^^e-^ "^ -2^yco,i3^ 
Jo 

/•oo 

dx / dyx"'-^y"'-^e-'''-y'-^''y''°"^. 
^0 



(B6) 



(B7) 



B.3. Smoothing with a top-hat filter 

We now consider a top-hat filter of radius unity and normalized 
to unity. Its Fourier transform is then 

W{k) = 3k-^{smk-kcosk) = 3, /yA;-^/^J3/2(A;), 

where Jp{x) denotes the spherical Bessel function of order p. 
We decompose equation (A23) in two terms according to 



(B8) where It h. can be computed analytically and Hi\i can not. That 

is 



We can now evaluate F-i(/3), by introducing the new variables 
V and w defined by 



dk dl 



k"l" 



k-^'^.hl2{k) 



X = V — w cot l3, y = 
The integral becomes 

F-i(l3)'- 



sin/3 



d(x,y) 



sin p Jo Jw cot/3 



d(v, w) 



dv e 



1 



sin/3 



(27r)*' " ' V" V 2 ^ 

l~^^^ -h/iil) \k + l\-^^^ .h/2(\k + l\) (flo + fl2 cos^ ' 

f k 
+a4 cos'* 9 + (ai cos 9 + a3 cos"' ^) ( T + 



H, 



th. 



dk dl 



k"l" 



In the polar coordinates v = p cos (p, w = p sin (p, the condition 
w > gives < (f < TT, while v > w cot /3 leads to cot (p > 
cot /3, oiQ < (f < /3. The integral is then given by 



(27r)*' ' ' V^V 2, 

k-"^hl2ik) l-y^ J3/2(0 \k + Z|-^/2 J3/2(|fc + l\) 



(1 



2 -.2X1 (k^ + fi) +X2klcost 



'■or 



k^ + P + 2klcosl 



F-i(l3)=\ 



P 



cos/3 



4 Vsin3/3 sin2/3^ 
Similar techniques can be used to evaluate 



By introducing the Legendre Polynomials Pm, and with the 
change of variable/ -^ —/justified further down, 7th becomes 



G-m = \ 



1 



/3cos/3 



/th^ = 



27 



De 

JG 



duk" 



4 Vsin2/3 sin3/3 

From the definitions (B7) and (B8), one can express Fn and 
Gn , for any value of n > — 1 , as 



i^„(/3) = (-2)- 



-n-1 



jn+l 



dcos"+i/3 



F-iiP), 



G„(/3) = (-2)- 



dcos"+i/3 



TT^G.m- 



Similarly, F-2 and G_2 are integrals over /3 of 2sin/3 times 
F-i and G_i respectively 

i^_2(/3)=-i/3cot/3, G_2(/3)=-^. 
2 2 sm /3 



167r3 V ^ JO JO J-i 

k-^'^.hl2{k) /-'/' J3/2(0 «;-'/' J3/2(W') 

(ao -Po(m) + ^2 -P2(m) + ctA Pa{u) 
k 

We have used 

M = cos(fc,Z), w = \k — l\, 

and 

02 04 3 

ao = 3(ao + -r + "t)) "1 = 6(ai H- -03), 

2 4 12 24 

"2 = 3(-a2 + -04), a3 = —03, a4 = ^704. 



(B9) 
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The term J3/2(w) is clearly a source of trouble because it mixes 
all three variables k, I and u. Fortunately one can use Gegen- 
bauer's addition theorem, which states for w = \Jk'^ + P — 2klu 
(e.g. Watson 1958); 



-''^J3/2(w) = V2^J2^-+m)(kir 



3/2 



m=0 



J3/2+m(k) =/3/2+m(0 —Pm+l(u). (BIO) 

Equation (B9) can now be expanded into an infinite series of 
products of integrals over k, I and u. Integral over u can be sim- 
plified using orthogonality properties of Legendre polynomials, 
to yield: 

oo ^1 , oo 

'^g(m) / duPq(u)—Pm+\(u) = 2'^g(2m + q). 

Let US define T{iy, 11,— \) = J^ J„{x) Jij_{x) x~^dx.ltvenfies 
(Gradshteyn & Ryzhik 6.574.2) 



T{v,^,-X) = 



r(A)r( "^^r^^' ) 



r^\ p/ —V+jJ,+X + \ \ p/ ^'+^+A + l \ p/ V — jJ,+X + \ \ 

for z/ + // + l>A>0. The second moment then is 

^ ' 47r 2 2 

and is finite for — 3 < n < 1. And the skewness is 



(BID 



/u 



{h^Y 



3 ao — 5 tti 



(1 - n)0 + n) 
(2 - n)(5 - n) 



with 

Tlniq) 
Qniq) 



1 

3&f 

+7^„(l)2(7(ao + a2)-9(al + a3)Qn(l)) 

OO 

+ ^Tln (qf ((4g + 3) (ao + a2 + 04) 

3=2 
-(4g + 5)(ai + a3)Qn(g))], 



J^(3/2,2g + 3/2,w-l) 
JP(3/2,3/2,n- 1) ' 
JP(3/2, 2q + 5/2,n)J^(3/2,2q + 5/2,n- 2) 



JP(3/2, 2g + 3/2, n - 1)2 
With the help of Eq. (B 11) and the identity 

r(x) ^ ^_^^,T(k + 1 - x) 



T(x-k) ' '' T(l-x) 
these ratios simplify to: 



TZn(q) = (-1) 



r(g + i±a)r(g + f)r(^)r(^) 



5-n^ 



Qn(g) 



r(^)r(f)r(g + ^)r(g + ^)' 

(1 - n)(2q + n)(2q + 3 + n) 
n(2q + 2 — n)(2q + 5 — n) 



One can now identify some very simple special cases: 



-3, 5; 



anal _ "0 

3 - U' 



1 canal 1. , 00-^1+02 «!. 

-2, ^3 =^(ao+ ^5 y), 



1 



n = 0, 5f"^' = ^(ao - ^). 






For n = 1, TZiiq) = (—1)*, Qiiq) = 0, and the series diverges. 
For other values of n the series has to be computed numerically 
and appears to converge very fast at least for n < 0. 

We now have to compute the extra term H in (i^"'). For 
n = —3, [S^) diverges and in the expansion of H according to 
Eq. (BIO) only the first term, with index m = 0, is divergent and 
contributes to skewness. One then gets ^3 = &J~^(ao + 8a;i/15). 
For all other values of n, H has to be computed numerically. 
This time, the integral cannot be reduced or simplified, and 
must be calculated as a three-dimensional one. Because of the 
presence of Bessel functions in the integrand, it oscillates very 
rapidly and we chose to use a Monte-Carlo integration. Ac- 
cording to the values of x\ and X2 (see Table) these integrals 
give a small contribution to the skewness and the numerical 
calculations could be done with a relative uncertainty that we 
estimate to be better than 1%. Thus, our numerical values of 
skewness are quite accurate, and the errors come mainly from 
the approximations we used for the logarithmic derivatives of 
the growth rates f\ and /2 when Vl ^ \. 

C. Measurement of ^3 : error-bars 

Here, we explain how we computed the error-bars shown in 
Fig. 5. Such errors are mainly due to finite volume effects and to 
the finite sampling of the measured density field. Note however 
that there will necessarily be some differences between various 
estimates of S^ia^) (corresponding to various snapshots in a 
simulation) at given n. One can estimate an average standard 
deviation associated to such variations 



Age 53 



1 



Nf 



ss 



1 ^^ 

snapshot S 



{S3,s - {s,)y 



(CD 



where #ss is the total number of snapshots considered (three in 
our case) and (^3) is the average of the measured S^ia^) over 
all the snapshots (53.5 is of course the measured value of ^3 in 
the snapshot S). Although this error is in fact correlated with 
the finite volume error and the finite sampling errors which we 
discuss hereafter, it can also be due to other defects, such as 
those related to initial conditions. We shall consider such an 
error as independent from others. 

To estimate the errors due to finite volume effects, we first 
fit the large-/9 tail of the measured distribution function P{p) by 
an exponential 



P(p) = -P(pmax) exp(-/3[p - Pmax]), 



(C2) 
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with 

/? = - 



1 



dP 



PiPmnx) dp 



(Pmax)- 



(C3) 



The error Apv due to finite volume effects on 53, is always in 
the same direction, i.e., 53 is always underestimated by direct 
measurements (CBS), since one misses the positive contribution 
due to the missing tail of the exponential. Following CBS, we 
use their Eq. (42) to write 



Afv5'3 



[(3 + 3^2^g^^g]^(Pi 



x) 



2«4 



0-2/3' 



(C4) 



with ( = iSpin^^. The fit is done using a least square fit method 
in the (p, log P) coordinate system (where an exponential looks 
like a straight line). 

The smoothing of the density field can be done either in 
real space or in Fourier space. The smoothing in Fourier space 
is much simpler to do than in real space, and much faster in 
terms of computer calculation. However, in the case one works 
in redshift space, it is impossible, practically, to have a sample 
obeying periodic boundary conditions, which makes the Fourier 
space calculation difficult, because of edge effects. Such arti- 
facts are expected to be small if the scale considered is small 
compared to the sample size, but it can become rather important 
(inducing errors larger than ~ 10%) as soon as £ ^ _Lbox/20. 
We thus decided to make all the calculations in real space, both 
for a top-hat smoothing and a Gaussian smoothing. In the case 
of a Gaussian smoothing, there is the added complication of the 
infinite spatial extent of the filter. One thus need to truncate the 
window at a finite radius r. We chose to make such a truncation 
at r = T£ with T = 3, i.e.. 



(C5) 



Wdr)xexp[-^] , r<3i, 



Wdr) = 0, 



> 31. 



A comparison of the results obtained using the procedure above 
with those obtained in Fourier space (without redshift mapping) 
shows very little differences, much less than one percent. More- 
over, we also checked that taking T = A instead of T = 3 did 
not significantly change the results. 

Note that the calculations are much more costly for a Gaus- 
sian filter than for a top-hat one, because, at a given site, we 
have to take into account the contribution of the density field 
up to 31 in the first case and only up to I in the second case. 
Hence, we did slightly less accurate calculations for a Gaussian 
smoothing than for a top-hat smoothing, in terms of number of 
sites sampled to evaluate the smoothed density distribution. In 
both cases, we computed the smoothed density field on a lattice 
with #s = ^g sites, but N^ was smaller for a Gaussian smooth- 
ing than for a top-hat smoothing. To sample appropriately the 
smoothed density field with this method, N^ should be large 
enough that the distance between two sites is smaller than the 
smoothing scale l considered. In other words, we should have 
I ^ Lbox/Ng, a condition fulfilled by our calculations. In order 



to evaluate the error associated with finite sampling, we used 
a different method for a Gaussian smoothing and for a top-hat 
smoothing. In the first case, we had an integral approach, and 
estimated the error due to finite sampling by assuming that we 
used the trapezoidal rule method to compute the integral. Let 
Sij^k be the estimate of the smoothed density contrast at the site 
(i, j, k). The estimator of (i^*^) is then 



E({6^)) = 



EC.- 



(C6) 



hj,k 



We estimate the error on this quantity by 

^FS {6'^) ~ ^ E [I^M-j,.l + \F2,,.k\ + \F3,,.k\] XC7) 



36N, 



hj,k 



with 



i.ij.k - "i+l,j,k '^"i-l,j,k 



■^2,i,j ,k 



"U+l,k ^ "i j-l,k 
^3,i,j,k = "ij^k+l "•" "i,j,k-l ~ ^"i,j,k 



26? 



(C8) 



In the case of a top-hat smoothing, the error was estimated 
in a more standard way, using the count-in-cell formalism. If 
we forget discreteness effects to simplify the notations (but we 
have fully corrected for them in the calculation of cP' and S^), 
we can write 



[Af5(^'^)f~-^[<^''^>-<^'^>'] 



(C9) 



This equation actually assumes implicitly that the cells were 
thrown at random in the sample, which is not the case in our 
algorithm. In other words, we certainly overestimate here the 



error, which should be proportional to N^ ' rather than N^ 



■1/2 



To estimate the error on S^ 
propagation of errors formula, i.e. 



((5^)/ ((^2) , we use the 



[Axx5'3]' 



[Axx(^^)f ^^ [(^^)Axx(^^)P 



with XX=FV or FS. After that, we finally write 

[A53]' = [AssS'3]' 



(CIO) 



(Cll) 



+ 77- E {[AFv5'3,5]' + [AFs5'3,5f}, (C12) 
snapshot S 

where [AxxS^s.s]^ stands for the estimate of the error of kind 
XX for the snapshot S. 
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